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ca ted 
ABSTRACT 
The purpose of this study was to investigate the effect 

of injectant stagnation temperature and molecular weight 
variation on the flow field generated from secondary injection 
of a gas normal to a supersonic stream. Experiments were 
conducted at a primary ateean Mach number of 2.80 in the 
Naval Postgraduate School supersonic wind tunnel. Experi- 
mental data have been correlated with various theories, 
showing both agreement and disagreement. Data presented 
include penetration of secondary jet into primary flow and 


shock shape as correlated with the second order blast wave 


theory. 
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CHAPTER I 
INTRODUCTION 


The interaction of a lateral jet with a supersonic stream 
is of fundamental importance to the understanding pf complex 
supersonic flows, The jet interaction application to thrust 
vector control, reaction control jets and fuel injection for 
supersonic combustion is of great engineering interest, 
Although extensive investigations!~ have been conducted of 
the supersonic gaseous injection phenomena, no systematic 
investigation of secondary jet stagnation temperature and 
molecular weight has been made. Most investigations have 
used sonic injection, so that in the absence of stagnation 
temperature and molecular weight variations it has not been 
possible to distinguish between the effects of jet momentum 
and jet mass addition, 

The objective of this thesis was be study the flow 
field induced by injecting a gas normal to a supersonic 
stream, where the primary interest was to discover the 
effect of injectant stagnation temperature and molecular 
weight in the flow field. These studies were performed in 
the Naval Postgraduate School supersonic wind tunnel at a 
Mach number of 2.80 by injecting gaseous nitrogen, argon 
and helium through a sonic nozzle mounted normal to a flat 
plate. The parameters varied during the experiments inelude 
injectant stagnation pressure and temperature, mass flow 


rate, molecular weight and injector exit area, By 


AA 





variation of injectant stagnation temperature and molecular 


weight, distinctions are made between jet mass and jet 


momentum effects. - 





CHAPTER II 
EXPERIMENTAL EQUIPMENT 


The injection equipment used in this thesis was arranged 
as in Figure 1. Compressed bottles of gaseous injectant were 
obtained commercially and fed into the system through a 
pressure regulator, The flow was then passed through a flow- 
meter, into the heater, and on to the sonic nozzle mounted in 
the flat plate. Photographs of the equipment assembly are 
shown in Figures 2 and 3. 

The method used to heat the secondary gas was simple, 


yet effective. Twenty-five feet of 3/8 inch 0.D. stainless 


os 
steel tube were coiled and placed in an insulated container 


as shown in Figure 4, Preliminary calculations indicated 
sonic velocity would be reached inside the heater tube at 
high temperature and pressure for standard tube diameters 
smaller than that used. An are welder was then connected 
across the ends of the coil, and the coil was heated by 
passing an electric current through the stainless steel tube. 
The injectant gas was subsequently heated through contact 
with the hot tube. The maximum injectant temperature was 
set at 1800°R due to the structural integrity of the heater 
at high pressures. | 

At the heater exit and just prior to the sonic nozzle 
entrance a device was installed (Figure 5) to decrease the 
velocity of the secondary flow to approximately stagnation 


conditions, The flow velocity was gently decreased through 
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a 5 degree diverging section until the cross-sectional area 
Was increased by a factor of ten, ensuring the velocity of 
the flow through the enlarged section to be less than 0.06 
Mach (or approximately stagnation conditions). At this 
point the temperature and pressure of the secondary fluid 
were measured, and these values were used in the isentropic 
equations for a sonic nozzle, A well rounded converging 
nozzle was then used to decrease the area ratio prior to 
entering the sonic nozzle. 

The sonic nozzle was a simple converging nozzle with 
parallel side walls at the exit section as shown in Figure 
6. Exit diameters of 0.1406 inches and 0.075 inches were 
used during the experiments. The sonic nozzle was calibrated « 
by means of an ASME sharp-edged orifice meter. The nozzle 
discharge coefficient was found to be 0.96 as summarized in 
the Appendix. 

A special phenolic supersonic tunnel block was fabri- 
cated to allow the flat plate pedestal to be mounted on the 
test section floor, The sonic nozzle could then be inserted 
through the bottom of the test section and through the 
pedestal without disturbing the flat plate. Thirty-nine 
static pressure taps were located on tne surface of the 
plate as sown in Figure 7. The generated pressures were 
measured with a mercury manometer. The side force data will 
not be discussed in this thesis, 


The NPGS supersonic wind tunnel is a blow-down model 


With a 4 inch by 4 inch test section, 6 inches long. The 





tunnel is capable of 5 minutes of continuous operation at a 
stagnation pressure of 50 psia. The Reynolds number at the 
injector port was 8.5 x 10° inches/inch for M = 2,8, ensuring 
a turbulent boundary layer in the region of the injector 
nozzle. 

The Schlieren system used had a standard off-set 
arrangement with collimating mirrors used in place of lenses. 
Tne mirrors, light source, and camera were positioned with 
the aid of a surveyer's transit, so that most of the parallax 
in the system was eliminated. With this arrangement, shadow- 
graphs could also be obtained by removing the knife-edge and 


changing the slit-source to a point-source, 


15 











CHAPTER III 
THEORETICAL ANALYSIS 


Two properties of the interaction phenomena, penetration 
height of the secondary fluid and radius of the induced bow ? 
shock, will be treated extensively in this thesis. 

The secondary gas is injected into a supersonic stream 
in the form of a highly underexpanded jet, as illustrated 
by Schetz and Billig? in Figure 8. The secondary jet pene- 
trates the primary stream to a certain height before it is 
completely bent over and proceeds downstream with the primary 
fluid. This penetration height depends on the properties of 
the two streams. Figures 9 and 10 illustrate interaction 
models by Zukoski and Spaida? and Hsia,° respectively. In 
each of these a penetration height is defined, and the e 
induced shock is noted. 

The radius of the induced bow shock can be treated as a 
function of the penetration height of the secondary fluid? 
or may be calculated from the second order blast wave analogy, 


as presented by Broadwell? and Dahm. 8 


Figure 11 shows 
Broadwell's flow model where only the radius of the induced 
shock is being considered. 

It is interesting to note that each of the flow models 
shown in Figures 9 through 11 was based on a blunt~body 
approximation of the interaction phenomena, but that the | 


results obtained by the model originators varied considerably, 


due to the basic assumptions made. 





The above models were used in this thesis as a guide 
for determining which parameters of the flow to vary and 


then were compared with the experimental results. 


PENETRATION HEIGHT 

Zukoski and Spaid,? in investigating the injection 
phenomena, obtained their scaling laws for secondary in- 
jection by setting up a model for a solid body, which would 
five a shock pattern similar to that produced by injection, 
By balancing the drag of the nose section of the blunt body 
against the momentum flux of the injectant, the character- 
istic dimension of the nose section could be found. They 
assumed the blunt body to be a quarter sphere and using 
Newtonian flow calculations found the force in the direction 


of the primary flow to be 


B= (r/8) wey,.¢.n?p,, (1) 


wo oe p 


where h is the radius of the quarter sphere and 


= ©, (Moo 9 Yeo ) 


is evaluated at the nose of the body. 

By injecting the secondary gas normal to a flat plate, 
the only contribution to the momentum in the direction along 
the plate will be that obtained from the expansion of the 
secondary fluid to the primary stream pressure. From the 


steady flow energy equation, 


hp, = hy + v,72 






or 


Vv, =2 °p, Ty), (l - T/T.) (2) 


Assuming expansion to the primary stream pressure P., 


the expression for the velocity vy can be written as 








2Y5 | al "is Yo : 2 
J ae ea. a 
or in terms of temperature, 
T 3 
2y; ° 
mj Ty -i od 9) P 1 +a Me (= 


The mass flow rate through a sonic nozzle with an 
effective exit area of CAy where C is the discharge co- 


efficient of the nozzle, is given by « 


* 1 4 g \ dati ae 
m4 = CA; Foil Ba vse) vi73 | (5) 


The force balance in the primary stream direction due to 





the momentum flux is then 
F =m, V (6) 


Equating (1) and (6), Zukoski and Spaid? solved for the radius 


h of the quarter sphere and found 





or in terms of stagnation conditions 





Mos A 
be. 2 Yo -t 2 \ély,, = 1) 04 Y} 2 
tre ay et z= Fe Py To “py 
2 2 Ya-L L — 8 
ai \yoar | ° hn = oe: er 
i J L+—s— Mm °3 


They then proposed that this calculated radius could be used 
as a measure of the scale of the disturbance produced by 
injection. Further, if the model is correct, it is expected 
that other flow features will also be proportional to changes 
in h, | 
feta” also assumed that the nature and extent of flow 

disturbances resulting from secondary injection are the same 
as those created by a blunt-nosed slender body. He then 
equated the nose drag of the slender body to the drag derived 
from momentum exchange between the primary and secondary 


streams and the effective force resulting from-energy ex- 


| change due to reaction and evaporation, giving 
ZA ° An 
ChA(A.V./2) = mV (12 =) 


where A is the projected area of the body in the direction 
OF flow, f,, and wv. are the density and velocity of the un- 
disturbed flow region, m, is the flow rate of the secondary’ 
fluid, and Ah is the energy exchanged per unit mass of the 
injectant due to chemical reaction and evaporation betweatt 


the primary and secondary fluids. Assuming that the shape 
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of the equivalent body was a quarter sphere followed by a 
half cylinder, Hsia found the nose radius to be given by 

r= (#-)"b+ 3) (a) 7 (9) 

mC), G et ie 
This can be thought of as the effective penetration height 
the secondary fluid would reach for a given flow rate of the 
secondary fluid. 
| Although both theories discussed above use a blunt-body 
model, the use of two different velocities leads to a fun- 
damental difference. Zukoski and spaid? assumed the secondary 
fluid expanded to the primary stream pressure, Po after 
leaving the injector nozzle, and therefore the final velocity 
attained was the velocity due to isentropic expansion to 
that pressure. Hsia® assumed the final velocity attained 
by the secondary fluid was that of the primary strean, Vv" 
This led te a penetration height by Hsia that was jet mass 
dependent and a penetration height by Zukoski and Spaid 
that was jet momentum dependent, 

The secondary mass flow from a sonic nozzle is pro- 
portional to the stagnation pressure, Fo, 3 and inversely 
proportional to the stagnation temperature, To , Of the 
secondary fluid. The momentum of the secondary fluid leaving 
the injector nozzle is only proportional to the stagnation 
pressure of the secondary fluid. The depth of penetration 
of the secondary fluid into the primary stream will be 
either jet mass or jet momentum dependent, or a combination 


3 


of both. Agreement between the Zukoski-Spaid~ and Hsia® 
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6 
penetration theory has been previously noted, but in the 
absence of stagnation temperature variation no distinction 


. can be made between the jet mass or the jet momentum effects. 


‘ BLAST WAVE ANALOGY 
The radius of the bow shock caused by secondary in- 
jection has been derived in terms of the second order blast 


8 and Broadwell? based their theories on 


wave equation. Dahm 
the blast wave analogy, where it was assumed that; 
(a) Injection of a secondary fluid causes the same pressure, 
velocity and density fields as a blunt-nosed body in 
a hypersonic strean, 


(b) The secondary fluid is accelerated to the primary 


stream velocity, V. 


' (c) The energy added per unit length of gas is equal to the 
change in momentum of the secondary fluid and the heat 
added to the primary fluid, 

(d) The blast wave analogy is valid throughout the entire 

region of the interaction of the primary and secondary 

streams, 

It can be noted again that gukoski and Spaid? also 
assumed a blunt-nosed body in a hypersonic stream (Newtonian 
calculations) to find the penetration height of the injectant, 

; but assumed the final velocity the secondary fluid reached 

was that from expanding to the primary stream pressure, Roy 

P instead of the primary stream velocity, vo 

The models of Broadwell’ and Dahm? were based on the 


blast wave analogy of Hayes.? According to Hayes' concept | 


: 





of hypersonic similitude, the pressure, density and velocity 
fields generated by an unyawed, axisymmetric body ina 
uniform hypersonic flow field are similar to the pressure, 
density and velocity fields generated by the explosion of 
a line charge. Once the pressure, density and velocity 
fields as a function of time t and distance from the line R 
have been obtained, tne corresponding fields for the 
axisymmetric body are derived by replacing t by x/V os where 
V,, is the velocity of the undisturbed primary flow and x 
is the distance along the axis of the flat plate from the 
apex of the detached shock, The concept of similitude is 
completed by relating the energy released by the line charge 
to the drag of the body. 

sakurai's?° second order solution for the shock wave 


radius R at any time t in the unsteady problem is given by j 


1 1 a s 
25 y F JoP. \ . 
R= £ . = ct ane 10 
Ltd oe? - sal 2E so co ( ) 
\ “Oe Lo ‘ 


where J Bw. hy are constants that depend on the specific 








heat ratio y,, (for y= 1.4, Jo = 0.88 and Ay = -1.989), PL 
and Cc. are the pressure and sound speed of the undisturbed 
primary fluid, and £— is the energy, measured per unit length, 


released in the creation of the shock wave. Assuming t = x/Vos 


and ct 


a, 1 
2 + a 2 
= y \a/ x \* Jo” = 


where x is the distance measured from the apex of the bow 


x/M,,» equation (10) reduces to 


shock along the centerline of the flat plate, V is tne 








velocity of the primary fluid at the injectant station, 


R= B/R,» xX = x/R, and 
3 
RS = @/enp,) . (12) 


is the characteristic radius related to the blast wave energy. 
As Broadwell’ pointed out, E is the energy per unit 
length in the completely symmetrical case, and for a flat 
plate all the energy 1s confined above the plate. Therefore, 
E/2 should be set equal to the total increase in energy 
due to injection. 
As stated in (c) the total energy added per unit length 
is equal to the change in momentum of the secondary fluid and 
the heat added to the primary fluid. Broadwell? derived the 


yey 


energy term as 


2 Ki 
" 2 ~1) M w 0 

B= 2m, Vi {1 + hE ae 2 2 a) (13) 
co 2 (y -21) uM. Ws % 


where the first term represents the energy related to the 
drag caused by the injectant fluids, and the second term 
represents the effect of the volume addition by the injectant 
gas. The development of the energy term by Dahm® parallels 
and extends that of Broadwell? but employs significantly 
different approaches in certain importane areas. Dahm 
derived the energy term for the second order blast wave 


radius as 
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2 T 
é 2+(Y -1)™ ° 
Yet, 1) H Ygra 2Y M- 5 7m 
2 2 | 
4 ¥, 7 M, [2 -1 
' (2 - 7) = oe cos | (4) 
Yoo j Os, 00 2+ {y,-1) M; 


where 6 is the angle between the surface of the plate and the 
axis of the injector (cos 6 = 0 for normal injection). The 
terms in the energy equation denote the increase in kinetic 
energy of the secondary fluid, a mass correction term, the 
heat added to the primary stream, and the effect of injecting 
at an angle. 

It is interesting to note that although Broadwell de- 
veloped the energy term, E, from an entirely different view- 
point than that of Dahm, his result agrees very closely for 
normal injection. 

Taking Dahm's value for E, for normal injection equation 


(12) reduces to 


aie a V3 24(¥-1)K i, an 
a a a a 
"Ye aici 0, ¢ * 27 Me es 


Expanding equation (15) in terms of the stagnation properties 
of the primary and secondary streams and assuming a constant 
primary stream molecular weight, specific neat ratio and Mach 


number, ti@h 











where A is a constant that depends on M . Since 


(17) 





further reduction may be made if for a particular secondary 
fluid botn the primary and secondary flow rates are held 


constant. Then, 


TS. 2 , 
i Rc |A+ BS (18) 
fo) T 
fe) 
> 0] 
: where B= oats) and A= A(M). 





ee LY aE A ae | 


CHAPTER IV 
EXPERIMENTAL RESULTS 


All parameters associated with the primary stream were 
held constant throughout these experiments, That is, air 
was used as the primary fluid, and the supersonic tunnel 
stagnation temperature and pressure, and tunnel Mach number 
were essentially cokatant for all runs. Therefore, when it 
is stated that the pressure ratio or mass ratio was varied, 
in effect, only the injectant pressure or mass was varied. 
The injectant stagnation pressure was varied from 12.5 psi 
to 300 psi, giving a range of pressure ratios from } to 6. 
The injectant stagnation temperature Was varied from room 
temperature (520° R) to 1590° R, giving a range of temperature 
ratios from 1 to 3. Nitrogen, argon, and helium were used 
as injectant gases to provide a substantial range for molec- 
ular weight effects and also to provide for any heat capacity 
ratio effects, The sonic nozzle exit diameter was also 
varied to change the mass flow rate of the injectant while 
keeping the preasivte and temperature constant. 

The individual experimental runs were conducted in the 
order shown in Figure 12. Schlieren photographs were taken 
of each run, and shadowgraph photographs were taken of many 
runs. In general the shadowgraph pictures taken did not add 
to the detail found in the Schlieren pictures and, therefore, 
were not extensively used. The static pressures displayed on 


the manometer board were also photographed at the same time 











the Schlieren photographs were taken. Polaroid type 55 P/N 
film was used for each Schlieren picture so that the negatives 
could be enlarged for comparison. The data taken from the 
photographs, which included penetration height of the in- 
jectant and radius of the bow wave, were obtained by either 
magnifying the photograph and measuring directly or placing 
the negative in an optical comparator. Seventy-seven 
individual runs were made for this thesis, each one varying 
the parameters of the injectant. Typical Schlieren photo- 
aegis obtained are shown in Figures 13 and 14. 

The experimental results are presented in two sections. 
Figures 15 through 20 depict the penetration heights obtained 
by varying each parameter of the secondary flow, and Figures 
21 through 31 illustrate the variation of the bow shock 


radius under the same conditions. 


PENETRATION HEIGHT 

The depth of penetration of the secondary fluid into the 
primary stream was scaled directly from the photographs and 
was taken to be the maximum height reached by the bent-over 
highly underexpanded jet. Figures 15, 16, and 17 show the 
penetration heights attained by argon, nitrogen, and helium 
respectively. In each figure the penetration height is 
plotted as a function of the stagnation pressure ratio, with 
the runs at various temperatures and nozzle exit diameters 
included, In each of these figures reference is made to the 
Zukoski and Spaid? theory for penetration height, showing 


agreement with their predictions, Although the stagnation 


2? 





temperature ratio was varied by a factor of three while 
holding the pressure ratio constant, no essential differences 
in tne penetration height could be detected. Increasing 
the stagnation temperature of the secondary fluid while 
holding the pressure constant causes the secondary mass 
flow rate to decrease. Since the penetration height of the 
secondary fluid was not affected by this action, either tne 
penetration height is a function of the secondary momentum 
or is at most a weak function of Moga! It should be 
noted that the slope of the line indicating the Zukoski- 
Spaid theory does not coincide exactly with that plotted from 
the experimental values obtained. This difference may be 
attributed to the exact value of the exponent of the pressure 
ratio term in equation (8). 

It may be noted in Figure 17 that considerably fewer 
data points are indicated, In general, it was found taat 
for the same quality Schlieren photograph the underexpanded 
helium jet was much more difficult to detect than was the 
argon or nitrogen jet. 

Figure 18 is a composite plot of Figures 15, 16 and 17, 
showing the mean penetration heights of the three gases as 
a function of the stagnation pressure ratio. As indicated 
there is a difference in the penetration heights of the 
three gases. Zukoski and Spaid's equation (8) is nota 
function of molecular weight; however, it does provide for 


heat capacity ratio variations, Argon and helium, naving the 


same value of heat capacity ratio, should have identical 





penetration heights under the conditions presented. Nitrogen, 
having a smaller value of heat capacity ratio, aheottia. have a 
penetration height 4 per cent higher than that of argon and 
helium, as calculated from equation (8). This is indeed the 
case when comparing the penetration heights of nitrogen and 
helium, but argon lies 10 per cent above nitrogen and 14 per 
cent above helium. This result shows that the penetration 
height is a function of the molecular weight. 


Equation (8) may be put in the form 


2Y < a rit cL P 
re ae ee per? 
h o (- | 5 nl iM ( > ) By <<l (19) 





Holding the primary stream parameters and the secondary mass 


flow rate constant, equation (19) can be further reduced to 


a 4 2 
h (5) (2, 15, | (20) 


Experiments were next conducted to investigate the 
penetration height while holding the secondary mass flow rate 
constant. For each gas considered the stagnation pressure 
and temperature were adjusted to give a mass flow rate of 
0.0065 lbs/sec ant 0.010 lbs/sec. The results are sosgentel 
in Figure 19. It is evident that for a given mass flow rate, 
the penetration height increases with an increasing stag- 
nation temperature ratio and also increases with a decreasing 
molecular weight. At a constant stagnation temperature of 
795°R, the following results were obtained from Figure 19, 


using equation (20): 





Calculated Measured 


h h 
No No 
hh. = 1.19 — = Tegaledl 
A n : 
A 
h nh 
— = 1.78 o = 1,89 ; 
A A 
nh h 
He He _ 
r— = 1-50 — = 1-70 
N2 No 


Disregarding the heat capacity ratio dependence of the 


penetration height, 1.e., the term 





produces: ‘ 
Calculated Measured } 
h h 
No N2 
—— 1.093 wo = ie 
A A 
h h 
He 
A "A 
h h 
— = 1.63 — = is70 
No Ng 


Disregarding the heat capacity ratio dependence of the 
penetration height can be viewed as changing the character- 
istic velocity of the secondary momentum. If the velocity at ; 
the sonic nozzle exit, v*, is used in equation (8) in place of 


Vy: then the results obtained from Figure 19, using equation 


(20), are: 








Calculated Measured 


— = 1.05 -—— = 1.11 


_He = le VO He st 1.89 
= 1.69 = 1.70 


The calculated values using v* are in better agreement with 


The data shown in Figure 19 scale as the 4 power of hy. 
When the heat capacity ratio dependence was neglected, the 
data also scaled approximately as the 4 power of the 
molecular weight, which is in agreement with the theory. 
The relatively large difference between the measured hops 
ratio and the expected value is attributed to a molecular 
weight effect. 

For a given gas (constant molecular weight) and the 
above stated conditions, 


2. 


oe ea 
"3 


There is an expected increase in penetration height of 32 
per cent between the temperature of 530°R and 1590°R for 
all gases.. Measured values found from Figure 19 were 44 
per cent increase for argon, 33 per cent for nitrogen, and 
12 per cent for helium. The lighter molecular weight gases 


show less variation with To, than predicted by theory. 


the experimental results than those using V5s 


SHOCK RADIUS 


Figures 20 through 22 show the shock radius variation 
for the various gases at a constant pressure ratio, varying 
temperature, The radius R was measured perpendicular to . 
the flat plate axis at distances x downstream from the in- 
jector axis. This presented a fixed reference point from 
which the actual shock radius could be measured. Comparing 
these figures, it is apparent that as the molecular weight 
of the injectant gas decreases, a greater variation occurs 
in the shock radius with temperature. 
Figure 23 shows the variation of the bow shock with 
molecular weight, while holding the stagnation pressure 
ratio and temperature ratio constant. Examination of 
equations (11) and (17) indicates that the radius of the : 
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bow shock should vary approximately as m *(a+ BR,)* for 


}es 


ee 


‘ 
these conditions, where A= A (Mm) and B= B vee Tt 


The variation of the bow wave with temperature and 
molecular weight was very evident when viewing the radius 


of the bow wave. When this radius is non-dimensionalized 





with respect to the characteristic radius, Ro» related to 
the blast wave energy, the effect is a convergence of all 
radii into a single line. This has been shown quite well 
in Figures 24, 25, and 26 respectively. The second order 
blast wave prediction, equation (11), has been added to 
these graphs to show the relationship between theory and 
measured values, In each case the measured values followed 


the pattern set by the second order theory with wnat appears 


, 


to be a constant displacement. However, when the mean lines 
from the three gases are placed on a composite plot, as in 
Figure 27, a molecular weight dependence is present. The 
spread between the three gases, as compared with any measure- 
ment error, is too large to be taken as experimental scatter. 
It may be noted that if a proper scaling radius, Ro based on 
E, is selected, then all curves should collapse to a single 
Line, 

Equation (18) was used as the basis for the next set of 
experiments, where the injectant mass flow rate was held 
constant and the radius R was viewed as a function of To ; 

Wy and Ys Figure 28 shows the variation of the shock radius 
with injectant stagnation temperature. Equation (18) was 
used to predict the increase of the radius with temperature 


with the following results: 





Calculated Measured 
R 
1590 ~ 1.10 =1590 = 1.13 
2530 2530 
BR R 
ee = 1,06 1080 = 1,06 
530 530 


Equation (18) was then used to predict the increase of the 
radius with molecular weight and heat capacity radio effects, 


as shown in Figure 29. The results were: 


Calculated Measured 
R RB 
He = 1,23 He 21.29 
R 
A A 
33 
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R R 

N N 

= = 1,04 221.06 
A Ry 

R 

_He _ 1.17 Bye Pg 
By ce 

No 2 


In general, equation (18) showed excellent agreenent for the 
constant secondary mass flow case. Figures 30 and 31 are 

the non-dimensional plot of Figures 28 and 29 respectively. 
Again the three separate radii collapse into a single line 
when put in the non-dimensional form of equation (11). Some 
aifferences due to molecular weight effects may be noticeable 


in Figure 3l. 
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CHAPTER V 
CONCLUSIONS 


At constant Poy? varying To » the penetration height 
theory of Zukoski and spaid? showed good agreement with 
experimental results. The penetration height also showed 
some dependence on molecular weight, which is contrary to 
theory. 

When the secondary mass flow rate was held constant, 
the penetration height varied with oy" This rules out any 
theory based strictly upon a jet mass dependence, as proposed 


6 in equation (9). 


by Hsia 

The penetration height scaled as (i) as predicted by 
all theories. 

Ignoring heat capacity ratio effects, the penetration 
height scaled roughly as (R,)4 which is in agreement with the 
jet momentum theory of Zukoski and Spaid. Assuming a 
different characteristic velocity of the secondary flow 
(using v" instead of Vy): good agreement was obtained with 
theory without having to ignore the heat capacity ratio 
effects (v*= (v2) for a sonic nozzle). If the heat 
capacity ratio effects were also ignored using this new 
Characteristic velocity .. then the results reduced to that 
obtained using te v" is therefore a better characteristic 
velocity for the flow process than V5- 

In general, the penetration height, as a function of 


T,., differed for different molecular weights. This is not 


° 
J 
in agreement with theory. When considered as a function of 
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Po » the experimental results for the penetration height 
sade good agreemnt with the theory of Zukoski and Spaid 
for both the constant Po, and constant mass flow runs. The 
conclusion is that the penetration height is not a function 
of T or is, at most, a weak function of Ty ° 

The radius of the blast wave measured from the experi- 
mental data was in close agreement with that predicted by 
Dahm. The blast wave theory was not able to predict the 
molecular weight variation accurately when the pressure and 
temperature of the system were varied (Figure 23) but showed 
good agreement when the secondary mass flow was held constant. 

The penetration height, h, was found to remain constant 
for varying temperature, and the bow shock radius, R, did not. 
Therefore, h would not be a good choice as a cnaracteristic 
dimension to represent the entire interaction phenomena. An 
attempt was made to non-dimensionalize the height of pene- 
tration with respect to the characteristic blast wave radius, 
Ro: with the results shown in Figures 32 through 34 Figure 
32 shows the variation of n/R, with respect to pressure for 
flows at various temperatures. The temperature spread is due 
to the fact that h does not vary with temperature, while a. 
does, Figure 33 indicates that as To ,/7o,, increases at 
constant mass flow, there is a divergence of the three gases, 
with respect to n/R This is because h does not vary as 
to increases while R. does vary as a function of temperature 
with another function of molecular weight included. Typical 


values are; 
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1 

2 
argon Rio [0.728 + 0.427 T/T, 

al 
nitrogen RX [0.728 + 0,848 T /70,)° 


au 
helium Ri « [0.728 + 4,28 T/T 


Figure 34 indicates the molecular weight effect present in 
the above terms at constant temperature. 
one statement that can be made of all three figures is 
that the slopes of their lines are near zero. This would 
indicate that a correction term in Dahm's energy equation 
might be enough to cause all the slopes to reach zero. 
Figure 35 shows the effect of non-dimensionalizing the 
bow shock radius with respect to the penetration height as 
formulated in equation (8). This can be compared with the 
a results obtained in Figure 26, which shows the bow shock 
non-dimensionalized with respect to the characteristic 
blast wave radius. The blast wave radius, Ro? appears to 
be a better characteristic dimension for the shock radius. 
In general, the comparisons between the characteristic 


dimensions, h and Ry» and the experimental data indicate 





that h would be a good engineering correlation parameter 
for secondary nitrogen and argon flows but would be less 
accurate for the lighter molecular weight gases such as 
heliun. 

Further, the data would support the contention that the 
heating of a lighter molecular weight gas would be more 


effective in increasing the area of interaction than the 
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heating of a heavier molecular weight gas. Since the aero- 
dynamic interaction force would be proportional to the 
interaction area, the heating of a lighter molecular weight 
gas would generate more side force than a heavier molecular 
weight gas. This advantage of increase in side force is 
particularly evident in the case of constant mass flow of a 
lighter molecular weight gas. 

Differences between the experimental results and the 
theoretical predictions due to molecular weight effects may 
be caused by the mixing processes of the secondary gas with 
the primary stream. 

In general, the concept of jet mass or jet momentum 


' dependence proved to be a useful tool in arranging the 


experiments to study the parameters involved, 
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APPENDIX I 


The sonic nozzle becomes a flow metering device whenever 
the ratio of the exit pressure to the isentropic stagnation 
pressure is less than the critical pressure, when this 


occurs, the mass flow rate through the sonic nozzle is given 








by 
aah 1 
* °o, Yq 2 ir iz 
e ”~ a Y - 
oa as WTo, BR, (a) J | (30) 


A stagnation temperature-pressure reservoir was installed 
just prior to the sonic nozzle which allowed measurement 
of the stagnation temperature and pressure used in the 
above equation. 

The sonic nozzle was calibrated by means of an ASME 
sharp-edged orifice flowmeter, which was built to ASME 
standards as set forth in Reference 11. The flowmeter was a 
Standard 2-inch diameter pipe with an orifice designed to 
give an orifice-to-pipe diameter retio g of 0.15. This value 
of g ensured that the Reynolds number through the orifice 
would always be greater than 10,000 for the range of mass 
flow rates considered. Vena contracta pressure taps were 
located across the orifice plate, and the static pressure and 
temperature within the flowmeter were measured at a point 6 
diameters upstream from the orifice. 

The pressure differential across the orifice was in the 
range of 0.5 to 1.5 psi for the mass flow rates considered 


during calibration, An attempt was made to measure this 
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differential with a pressure transducer, but none were 
available that would give the precision required at these 
low differentials. <A mercury U-tube manometer was then 
fabricated by milling % inch semicircular slots in 1 inch 
thick plexiglass and bonding the two halves together. The 
result was a 2~inch thick plexiglass U-tube manometer 
capable of withstanding 300 psi pressure, with no apparent 
effects from creep. This performed satisfactorily during 
calibration, and pressure differentials could be read to 
a precision of 0.05 psi. The plexiglass U-tube manometer 
Was calibrated against a standard glass-tube manometer, 
with no differences detected. 

The mass flow rate through a sharp-edged orifice flow- + 


meter is given by 
| : 
— ? 
m= YKA, [2 /, ap | 


where Y (Expansion Factor) is determined by 


2 , 
Fed wth 1295 (A,/A, ) Nts /Y4P, (31) 
for vena contracta taps, 


K = flow coefficient 
= 0.6133 to 0.5995 for Re = 10,000 to 100,000 
and g = 0.15 
A, = Orifice area 
Ay = pipe area 
Py = density upstream of the orifice 
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SP' = pressure differential across the orifice 


as) 
I 


pressure upstream of the orifice 


Table I shows the values obtained during calibration of the 
0.1406 inch diameter sonic nozzle. 

The sonic nozzle discharge coefficient, c, was defined 
as the ratio of the actual nozzle flow (flow measured through 
the sharp-edged orifice) to the flow calculated from the 
isentropic laws for a Sonic nozzle, 

= aot 
ISENTROPIC FLOW RATE 

For well-designed nozzles with straight axes having 
"pipe" Reynolds numbers measured at the minimum area of 10° 
or more, the discharge coefficient is of the order of 0.99.17 
The pipe Reynolds numbers during calibration ranged from 


+ to 1.4 x 105, and therefore slightly less discharge 


5:4 x 10 
coefficients were obtained, 

All nozzles used were checked to determine that their 
discharge coefficients fell in the same range as the one 


calibrated. 
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Figure 11. Broadwell Flow Model. 
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i Figure 3C. Non-dimensional Flot of Figure 28. 
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Figure 32. Conparison of Characteristic Disensions--Varying Tewperature. 
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Pigure 33. Comparison of Cnaracteristic Dinensions--Mass Flow Constant. 
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